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Okayama  University,  Tsitshimanaka  3-TI,  Okayama  700-8530,  Japan 


Abstract.  The  structures  of  Yukawa  systems  spherically  confined  under  the  condition  of  micro- 
gravity are  analyzed  by  molecular  dynamics  simulations.  It  is  shown  that  the  system  is  separated 
into  the  outer  and  the  inner  parts  and  well-developed  shells  appear  in  the  outer  part  at  low  tem- 
peratures. In  the  case  of  Yukawa  mixtures,  there  occurs  the  separation  of  species  according  to  the 
magnitude  of  charges. 


INTRODUCTION 


Dusty  plasmas  have  provided  us  with  an  important  example  of  nonneutral  plasmas.  One 
of  main  differences  from  other  nonneutral  plasmas  composed  of  ions  or  electrons  is  that 
dust  particles  have  large  masses  and  the  effect  of  gravity  plays  a key  role  in  determining 
their  behavior. 

Recently,  some  experiments  under  the  condition  of  the  microgravity  have  been  per- 
formed [I].  In  this  case,  the  gravity  can  be  mostly  neglected  and  we  may  observe  var- 
ious new  aspects  of  dusty  plasmas  which  have  been  hidden  by  the  effect  of  gravity  in 
experiments  on  earth.  We  may  also  expect  the  investigation  of  dusty  plasmas  under  mi- 
crogravity might  lead  to  some  new  method  to  synthesize  and  process  useful  materials. 
The  purpose  of  this  paper  is  to  extend  our  theoretical  and  numerical  analyses  on  the 
structures  of  dusty  plasmas  under  the  gravity [2][.l]  and  analyze  dusty  plasmas  under 
microgravity. 

We  assume  that  the  particles  interact  via  the  isotropic  repulsive  Yukawa  potential 


where  q is  the  charge  and  1/k  is  the  screening  length.  It  is  known  that  the  ion  flow  in 
the  sheath  or  other  effects  may  induce  ani.sotropy  into  the  interaction.  Here  we  show 
interesting  features  of  the  system  without  those  effects,  expecting  that  they  remain  so 
when  the  anisotropy  is  taken  into  account  as  a perturbation. 

In  the  case  of  uniform  one-component  Yukawa  system,  we  have  two  parameters, 
the  density  n and  the  temperature  T,  which  we  characterize  by  two  dimensionless 
parameters 


r = 


kgT  a 


and  ^ = Kfi, 


(2) 


where  n is  the  mean  distance  defined  by  a = 
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In  the  case  of  mixtures  with  species  / = 1.2....,  we  have  the  densities  the  charges 

<1,,  the  masses  ni,,  and  the  total  number  density  n,„i  as  parameters:  n,-/  = "i-  We  here 

consider  the  case  of  binary  mixtures  and  denote  the  species  by  i = 1,2.  We  delinc  the 
parameter  F i (denoted  as  F in  figures)  by 


r,  = 


<ii 

l-nTa' 


O) 


where  the  mean  di.stance  n is  defined  for  the  total  density  n,„,  as  a = (3/-l7rn,„,)''^*.  In 
thermal  equilibrium,  static  properties  arc  characterized  by  the  following  four  parameters 


r,- 


l>2 

"I 


(4) 


EXTERNAL  CONFINING  POTENTIAL 

We  consider  Yukawa  .systems  in  the  environment  of  very  small  gravity  and  assume  that 
Yukawa  particles  are  spherically  confined  by  an  external  potential.  As  the  source  of  the 
external  potential,  we  assume  that  the  sheath  with  the  uniform  space  charge 

c/i.o,  (/•>»■(,)  (5) 

surrounds  the  domain  of  electrically  neutral  plasma  of  radius  7o.  We  denote  the  clee- 
trostatic  potential  due  to  the  charges  in  the  sheath  by  <;'),.,/(/■).  The  external  potential  for 
Yukawa  particles  with  the  charge  -</  is  given  by  or 

= ('■<  '■())  (6) 

and 

, , , r/(4/rr/).s/,ri!)/3 

-'/W  .,/('■)  = 

'II 


+ 


('■  > 'o)- 


(7) 


NUMERICAL  METHOD 

In  order  to  have  as  many  particles  as  possible  in  the  system,  we  here  adopt  the  molecular 
dynamics  with  the  (){N)  .scaling  (3,  6J.  The  typical  system  size  in  our  simulations  is 
N = 10'.  The  temperature  is  controlled  by  the  Nose-Hoover  method.  For  mixtures, 
we  assume  the  same  mass  for  all  species.  In  the  classical  case,  the  static  properties 
do  not  depend  on  particles  masses,  though  the  relaxation  time  in  general,  especially  for 
the  energy,  increases  with  the  mass  ratio:  In  this  paper,  we  are  interested  in  only  static 
properties  after  relaxation. 
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RESULTS  AND  DISCUSSION 
One-Component  Yukawa  System 

Some  examples  of  the  results  of  simulations  with  10'^  particles  are  shown  in  Figs.  1, 2, 
and  3,  where  the  radial  distributions  and  the  position  of  shells  are  shown.  The  results  are 
summarized  as  follows.  ( 1 ) The  system  is  separated  into  the  outer  and  the  inner  parts.  (2) 
In  the  outer  part,  a thin  shell  is  formed  at  the  surface  of  the  system.  The  distribution  of 
particles  on  this  shell  is  a two-dimensional  lattice  with  defects.  Inside  of  the  outermost 
shell,  there  appears  a series  of  shells  as  in  the  case  of  one-component  plasma  confined 
by  parabolic  potential.  The  average  density  is  higher  than  the  inner  part.  (3)  The  shells 
become  higher  and  sharper  with  the  decrease  of  the  temperature  for  fixed  values  of  At 
fixed  temperature,  they  become  higher  and  sharper  with  the  decrease  of  ^ or  the  increase 
of  the  coefficient  of  confining  potential.  (4)  In  the  inner  part,  the  average  density  is  low 
and  the  shell  structure  is  not  clear. 


Yukawa  Mixtures 


The  results  for  the  mixture  of  Yukawa  particles  where  the  total  number  is  1(F  and 
rhilui  = 1,  q^i/qx  = 10,  and  q-z/q\  = 5 are  shown  in  Figs.  4,  5 and  6.  One  of  the 
most  remarkable  phenomena  may  be  the  separation  of  particles  with  different  values 
of  charges.  This  tendency  becomes  strong  with  the  increase  of  the  screening  length  and 
the  decrease  of  the  temperature.  As  for  the  relative  position  of  shells  of  each  species,  we 
are  still  trying  to  give  simple  interpretations. 


Analysis  of  Structures  at  Low  Temperatures  in  Fluid  Limit 


Here  we  present  a theoretical  analysis  on  the  structure  of  one-component  Yukawa 
particles  at  low  temperatures  based  on  the  variational  approach  which  is  successful  for 
the  two-dimensional  case  [7].  In  what  follows  we  adopt  the  fluid  approximation  and 
express  the  density  of  particles  by  p(r).  The  total  number  of  particles  N is  related  to 
/;(r)  by 

J drp{r)  = N.  (8) 

The  interaction  is  expressed  as 


// 


(Irdt 


,exp(- 


r-r' 


p{r)p{r') 


(9) 


and  the  external  potential  is  given  by 

= (-7)  |dr0er((i’)p(r). 


(10) 
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FIGURE  1.  Distribution  vs.  radial  distance 

measured  in  a in  spherically  confined  Yukavna  FIGURE  4.  Distribution  vs.  radial  distance 

system  under  microgravity:  Effect  of  screening  measured  in  a in  Yukawa  mixture  with  charge 

length  for  r = -50.  ratio  5 (P  = Ei). 


FIGURE  2.  The  same  as  Fig.  1:  Effect  of 
temperature  for  ^ = 0.1. 


0.10  I 1.0 


FIGURE  3.  Position  of  shells  measured  in  a 
vs.  f . 


FIGURES.  The  same  as  Fig.  4 for  ^ = 1.0. 


FIGURE  6.  The  same  as  Fig.  4 for  charge 
ratio  10. 
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Let  us  find  p{r)  which  minimizes  the  value  of  Ui„t  + U„i  under  the  condition  (8).  The 
variation  with  respect  to  p(r)  leads  to 


f dr  {-q)(j)ext{r)  + q‘^  1 

(11) 

J [ J 

and 

j 

|r  — r 1 

drdp{r)  = 0. 

(12) 

Denoting  the  Lagrange’s  multiplier  by  p,  we  have 

(-g)</>aTt(r)  + <7^j 

(13) 

and  the  value  of  //,  is  determined  by  (8). 

When  we  denote  the  source  charge  density  of  the  external  potential  by  ep^xt{r)  as 

Ai^pxt(r)  = -47rcp„,(r),  (14) 

we  have  from  (13) 

2 -2 

When  the  external  charge  density  is  given  by  (5),  we  have 

^2 

= (»'<»•())  (16) 
and 

We  observe  that  the  density  is  low  in  the  central  part  r < r<)  and  increases  stepwise  by 
(c/g)ns/,  at  r = tq,  separating  the  inner  and  outer  parts  of  the  system.  This  tendency  is 
in  agreement  with  the  behavior  of  the  results  of  simulation  when  the  layering  structure 
in  the  outer  part  is  averaged;  The  latter  seems  to  appear  when  the  effect  of  discreteness 
or  the  correlation  is  taken  into  account. 

When  we  denote  the  maximum  radius  of  distribution  by  r,„,  the  condition  (8)  gives 


This  relation  and  the  condition  p(r,„)  = 0 determine  the  parameters  /<  and  I'm. 
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Discussion 


The  average  distribution  may  be  understood  within  the  consideration  in  the  fluid  limit 
and,  since  the  repulsion  between  particles  increases  with  the  charge,  we  may  naturally 
have  the  separation  of  charges  in  the  case  of  mixtures.  The  formation  of  structures 
beyond  the  average  distribution  is  related  to  the  correlation  between  particles  [2,  3]. 

It  has  been  shown  that  particles  distributed  in  the  inner  part  may  serve  as  the  source 
of  the  force  in  the  radial  direction  which  plays  the  role  of  the  gravity  in  the  formation 
of  layered  structures  [4].  The  shells  on  the  periphery  in  the  one-component  case  thus 
correspond  to  the  one-dimensional  layers  under  the  gravity.  The  structures  in  the  case  of 
mixtures  seem  to  need  more  elaborate  analysis. 


CONCLUSIONS 

In  this  paper,  we  have  analyzed  the  behavior  of  particles  interacting  via  the  repulsive 
isotropic  Yukawa  potential  when  the  effect  of  gravity  is  neglected.  In  the  confined 
Yukawa  system  of  one  species,  the  layered  (sliell)  siructure  is  formed  from  outer  bound- 
ary of  the  .system.  This  structure  becomes  clear  and  shatp  with  the  deerease  of  the  tem- 
perature or  the  increase  of  the  screening  length.  Within  each  layer,  particles  are  orga- 
nized into  two-dimensional  lattice  with  defects. 

When  we  have  two  kinds  of  Yukawa  particles  in  our  confined  system,  we  have  a 
separation  according  to  the  relative  magnitude  of  charges:  The  species  with  larger  charge 
tends  to  distribute  on  the  periphery.  There  is  also  the  tendency  to  form  layered  structures. 
Both  becomes  significant  when  the  temperature  is  low  and  the  screening  length  is  large. 
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